Enhancement of ferroelectricity at metal/oxide interfaces 
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By performing first-principles calculations on four capacitor structures based on BaTiOs and 
PbTiOs, we determine the intrinsic interfacial effects that are responsible for the destabilization 
of the polar state in thin-film ferroelectric devices. We show that, contrary to the established 
interpretation, the interfacial capacitance depends crucially on the local chemical environment of 
the interface through the force constants of the metal-oxide bonds, and is not necessarily related to 
the bulk screening properties of the electrode material. In particular, in the case of interfaces between 
AO-terminated perovskites and simple metals, we demonstrate a novel mechanism of "interfacial 
ferroelectricity" that produces an overall enhancement of the ferroelectric instability of the film, 
rather than a suppression as is usually assumed. The resulting "negative dead layer" suggests a 
route to thin-film ferroelectric devices that are free of deleterious size effects. 

PACS numbers; 71.15.-m, 77.84.-s, 73.61.Ng 
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Capacitors based on ferroelectric perovskites are po- 
tentially attractive for applications in nanoelectronics, 
such as non-volatile random-access memories and high- 
permittivity gate dielectrics. Thin-film geometries are 
sought after for optimal efficiency and information stor- 
age density [1]. However, in the thin-film regime, strong 
size-dependent effects arise. These can considerably 
worsen the attractive functionalities of the ferroelectric 
material by, e.g., reducing the dielectric response [1, Q 
and causing rapid polarization relaxation Under- 
standing and addressing these deleterious effects is cru- 
cial for future progress [5|. 

Several hypotheses have been formulated to interpret 
size effects in thin-film ferroelectric capacitors. Usually 
the experimental data can accurately be described in 
terms of the "series capacitor model," where the over- 
all capacitance density C of the device is written as 
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Here Cb is the bulk permittivity of the dielectric material 
of thickness t and Ci_2 are capacitances associated with 
low-permittivity layers located at the film/electrode in- 
terfaces. Such interfacial capacitances produce a depo- 
larizing field 0, 0] that strongly reduces the dielectric re- 
sponse in the paraelectric regime and can destabilize the 
single-domain ferroelectric state. An accurate knowledge 
of the interfacial capacitance and its dependence on the 
combination of ferroelectric and electrode materials is, 
therefore, crucial for device design. 

Interestingly, such interfacial "dead layers" are present 
even in high-quality epitaxial systems where the de- 
fect density is very low This suggests possible 
fundamental issues affecting the polarization at the 
metal/ferroelectric boundary. In particular, within the 
Thomas-Fermi model, many authors focus on the imper- 
fect compensation of the polarization charges due to the 
finite electronic screening length of metallic electrodes. 



To explain the substantially better performance of SRO 
electrodes as compared to Pt ones, as has been observed 
experimentally ^, i3(| , one can then invoke the lattice con- 
tribution to the screening, which is expected to be sig- 
nificant 0] in oxide electrodes. 

The above interpretation, however, is clearly unsat- 
isfactory at the microscopic level. The region where the 
interface effects occur is as thin as a few interatomic spac- 
ings . In this region the local chemical and electrostatic 
environment departs significantly from that of either par- 
ent material, and a description of the interface in terms 
of bulk parameters is unjustified [1]. To describe such 
effects, a full quantum-mechanical treatment is required. 

First-principles calculations have already been invalu- 
able in understanding the properties of nanoscale ferro- 
electrics. Several studies have focused on the effects of 
strain [l^, HI and interfacial electrostatics [1, [H, 
in determining the spontaneous polarization of symmet- 
ric [14)] or asymmetric [15.] capacitor geometries with vari- 
ous ferroelectric/electrode combinations. While most au- 
thors agree on the existence of a depolarizing field that 
tends to suppress single-domain ferroelectricity, a clear 
description of the microscopic mechanisms causing it has 
not yet emerged. Until it does, it will be difficult to over- 
come this deleterious behavior. The recent development 
of rigorous first-principles methods to define and con- 
trol t he p olarization in a metal/ferroelectric heterostruc- 
ture [l6l [13] now provides the opportunity to perform 
such analysis. 

In the present work we concentrate on ultrathin ferro- 
electric capacitor structures in which the ferroelectric is 
PbTiOa (PTO) or BaTiOa (BTO) and the metallic elec- 
trode can be SrRuOa (SRO) or Pt. Aside from being 
technologically relevant combinations, these are diverse 
enough to cover several kinds of behavior while remaining 
focused enough to allow an in-depth analysis. We find 
that the interfacial dielectric response is more complex 
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than usually assumed in phenomenological models, where 
only the penetration of the electric field in the electrode 
is considered. In particular, we develop a rigorous theory, 
based on the modern theory of polarization, in which the 
interface-specific effects of purely electronic screening and 
of interatomic force constants are both taken fully into 
account in assessing the overall performance of the capac- 
itor. Based on our analysis, we demonstrate a covalent 
bonding mechanism that yields a ferroelectric behavior 
of the interface between AO-terminated films and simple 
metals. In this case one finds an overall enhancement 
of the driving force of the film towards a polar state, 
rather than a suppression as is predicted by semiclassical 
theories. 

Our calculations are performed within the local- 
density approximation of density-functional theory and 
the projector-augmented-wave method [l8| as imple- 
mented in an in-house code. An isolated capacitor with 
semi-infinite leads is modeled by a periodic array of al- 
ternating metallic and insulating layers, where the thick- 
nesses of the metallic slab is treated as a convergence 
parameter and the ferroelectric film is N unit cells thick 
(the actual values are reported in Table |T| . We focus on 
symmetric interface terminations of the SrO/Ti02 type 
for SRO electrodes, and of the Pt2/A0 type for Pt elec- 
trodes. In all cases we constrain the in-plane lattice pa- 
rameter to the theoretical bulk SrTiOs value. We use 
our recently developed fixed-Z? [l7[ method, combined 
with the extensions to metal/insulator heterostructures 
of Refs. and ITqI . to analyze the linear response of 
the paraelectric structure to a small polar perturbation. 
This regime is relevant for the existence of a (meta) stable 
single-domain ferroelectric state, and allows for a micro- 
scopic analysis of the physical ingredients contributing to 
the screening, without complications arising from asym- 
metries or pathological band alignments . 

In order to assess the stability of a given structure 
against a polar distortion, we are concerned with C~^, 
its inverse capacitance per in-plane unit cell, which can 
be expressed as [13] 
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where U is the internal energy per in-plane cell, S is the 
cell area, and V is the potential drop across the capacitor 
plates. We note that is a "generalized" inverse ca- 
pacitance that is meaningful even when it is negative [13] , 
in which case it signals the appearance of an instability to 
a ferroelectric distortion in the capacitor arrangement - 
the more negative, the stronger the driving force towards 
a polar state. 

We calculated for the four capacitor heterostruc- 
tures considered in this work using a finite-difference ap- 
proach. In particular, for each system we performed two 
complete electronic and structural relaxations, first at 
-D = by imposing a mirror symmetry, and then at 
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(mVF) 


iVcrit 


Aeff (A) 


BTO/SRO 


6.5 


-1.553 




2.280 


4.85 


0.202 


PTO/SRO 


6.5 


0.439 




1.727 


7.45 


0.153 


STO/SRO 


6.5 


3.876 




1.647 




0.146 


PTO/Pt 


8.5 


-1.427 




1.258 


5.42 


0.111 


BTO/Pt 


8.5 


-7.920 




0.037 


0.08 


0.003 



TABLE I: Calculated inverse capacitance densities C~^S, in- 
terfacial inverse capacitance densities Cf^S, critical thickness 
A'^crit and effective screening length Aefi for the four capacitor 
heterostructures discussed in this work. Data for SRO/STO 
from Ref. [l^ are reported for comparison. 



D ~ 0.001 a.u. (a value that is small enough to ensure 
linearity) using our finite-field technique 
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possible because the fixed- Z? method allows one to define 
and compute the equilibrium state of an insulating sys- 
tem, within a given set of symmetry constraints and at a 
given value of macroscopic D, even in configurations that 
would be unstable at fixed electric field (e.g. a ferroelec- 
tric near the saddle point of its double- well potential). 
Finally, we extracted the induced bias potential and dis- 
cretized Eq. ([2]) to obtain the inverse capacitances, whose 
values are reported in Table [H The positive sign of 
in the PTO/SRO case indicates that the system is para- 
electric for the considered PTO thickness, while the other 
structures are in the ferroelectric regime. 

A quantity that has received considerable attention in 
the recent past [T^ . 13 . 14 1 is the critical thickness 



for ferroelectricity, defined as the minimum thickness for 
which a polar state exists. The usual strategy to deter- 
mine this quantity computationally is to perform several 
calculations for varying thicknesses until a spontaneous 
polarization appears. Instead, we find it much more con- 
venient to exploit the power of the fixed- Z? approach and 
obtain such information directly from a calculation on a 
single thickness. 

To that end, it is useful to look first at the spa- 
tial decomposition of the inverse capacitances shown in 
Fig. [l] where we plot the local inverse permittivity pro- 
files e'^{x) = d£{x)/dD (evaluated at D=0), where £ is 
the X component of the y-z-averaged electric field along 
the stacking direction x. (These are calculated as in 
Ref. [13, except that our use here of the fixed- Z? method 
allows treatment of cases having < 0.) The profiles 
are generally characterized by two regions, an interfacial 
part where the imperfect screening manifests itself as a 
positive peak in e~^{x), and the deep interior of the in- 
sulating film where e^^{x) converges to the bulk value. 

This naturally suggests a local decomposition between 
interface and bulk effects, where the overall stability of 
the centrosymmetric state emerges from the competi- 
tion between these two usually opposite contributions. 



Rewriting Eq. ^ as C"! = Cf ^ 
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FIG. 1: Inverse local permittivity profiles for symmetric ca- 
pacitor structures studied in this work. (For the SRO elec- 
trode, the STO curve from Ref. is added for comparison.) 
Inverse permittivities calculated for bulk STO, BTO, and 
PTO are shown as dashed horizontal lines. Electrodes and 
oxide layers closest to the interface are indicated by shading. 
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Pt/PTO PTO /SRO SRO /BTO BTO/Pt 




FIG. 2: Plots of Callen dynamical charges, arranged head-to- 
head and tail-to-tail to emphasize rapid convergence to com- 
mon bulk values. Colored symbols denote A-site cations (red 
squares), B-site cations (light green diamonds), AO oxygens 
(Ol, blue down triangles), and BO2 oxygens (02, dark green 
up triangles). The Pt atoms (black symbols) are labeled ac- 
cording to their projection onto the interface plane. Of the 
four independent Pt sites, one is shared with the A cations 
(squares), another is shared with Ol (down triangles) and the 
remaining two, equivalent by symmetry, are shared with 02 
(circles). 



obtain for a symmetric capacitor 
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(3) 



Here N is the number of bulk cells and C^^ — inch/iebS) 
is the bulk inverse capacitance per unit cell (cb is the cell 
parameter along the field direction). We take Eq. ([3]), 
based on a symmetric capacitor, to be our definition of 
the interfacial inverse capacitance. (Note that the defini- 
tion depends on the precise convention for specifying N; 
we adopt the convention that N is the nominal thickness 
as illustrated by the examples in Table HI) Then, given 
and C^^, we again use Eq. ([1]) to predict the critical 
thickness 
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for ferroelectricity in a (possibly asymmetric) capacitor, 
this being the value of N that yields an overall vanishing 
inverse capacitance. For an interface to a ferroelectric 
material {C^^ < 0), it is also natural to define the effec- 
tive dead-layer thickness Ndcad,i = ^C^^ /C^^ , in terms 
of which 



-/Vcrit — -/Vdcad.l + A^dcad,2 



(5) 



or, for a symmetric capacitor, iVcrit — '^Ndca.d- 

The calculated values for and A'^crit are reported 
in Table HI together with the values already calculated 
in Ref. [ig for paraelectric SrTiOa (STO) between SRO 
electrodes. Since many authors discuss interfacial effects 
in terms of effective screening lengths Aoff = C^^S/An 



rather than capacitances, we also report these values in 
the table. 

It is rather surprising that the calculated iVcrit of 
BTO/SRO is smaller than that of PTO/SRO. Indeed, 
PTO has a much larger spontaneous polarization Ps in 
the tetragonal ground state than BTO, and one would be 
tempted to think that the former is a "stronger" ferro- 
electric and therefore should have a smaller critical thick- 
ness. The calculated trend in iVcrit is even more surpris- 
ing on observing that of PTO/SRO is smaller than 
that of BTO/SRO, which would favor the opposite trend. 
It is clear from Eq. ([4]), however, that |Cj7^|, rather than 
Ps , is the parameter that defines the strength of the bulk 
ferroelectric instability; at the fixed STO in-plane lattice 
constant, our calculated \C^^\ is about twice as large in 
BTO than in PTO, which agrees nicely with the trend in 
A'crit- This highlights the fact that |C,^^|, rather than Pg, 
is the important parameter to consider when discussing 
the competition of bulk and interfacial electrostatics in 
ferroelectric capacitors. 

In the Pt-based systems considered here, AeS is sys- 
tematically smaller than in the SRO-based systems, in 
stark contrast with the predictions of the semiclassical 
models of Refs. [l^. (First-principles calculations have 
already demonstrated, in the case of STO paraelectric 
capacitors, that Pt electrodes can be intrinsically supe- 
rior to SRO electrodes Here, while in PTO/Pt 
Aoff is still significant, in the case of BTO/Pt it practi- 
cally vanishes, or equivalently, the dead-layer thickness 
of iVdead = 0.04 layers is almost zero! This suggests that 
novel effects must take place at this interface that go be- 
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System 


d(A) 


Spt (A) 


<5ao (A) 




BTO/Pt 


2.26 


0.46 


0.12 


O 


PTO/Pt 


2.11 


0.34 


0.19 





Unsupported 



TABLE II: Computed structural parameters for Pt2-AO in- 
terfaces. See Fig. [3]for the definition of the symbols. 



yond the usual Thomas- Fermi arguments, which always 
predict the same, positive value independent of the fer- 
roelectric material. 

To investigate the origin of such effects, we perform a 
microscopic analysis of the electrostatic response to a I? 
field in the two Pt-based capacitors. The static inverse 
capacitance can be decomposed as 



c- 



(6) 



where {C°°)~^ is the purely electronic, frozen-ion value 
and —A is the lattice contribution. The latter can be 
expressed in terms of the Callen (or longitudinal) dy- 
namical charge tensor 



constant matrix K^, as 



2l| , and the longitudinal force 



V s 



1 r^L 



(7) 



where x is again the stacking direction. Note that the 
sums in Eq. ([7]) run over all atoms in principle, includ- 
ing those located deep in the metallic electrodes, but in 
practice these do not contribute because their dynam- 
ical charges are zero (see Fig. [5]). The convenience of 
using Eq. ([7]) instead of the usual decomposition involv- 
ing transverse quantities, e.g. Eq. (8) of Ref. [g^, is that, 
in a layered heterostructure, is short-ranged in real 
space [23|,|24|. This means that local bonding effects can 
be unambiguously separated from the long-range electro- 
static interactions, thus allowing for a detailed analysis 
of the microscopic mechanisms contributing to the polar- 
ization. 

It is clear from Eq. ([7]) that two basic ingredients con- 
tribute to the ionic polarizability, the force constants con- 
tained in K^, and the dynamical charge associated with 
a given degree of freedom, . Since the inverse of 
enters the sum, it is obvious that less stable bonds will 
yield an enhanced response, provided that the participat- 
ing atoms carry a significant Z^. It is reasonable then to 
expect that the relative weakness of the bonds formed by 
the ferroelectric films with Pt, compared to those formed 
with the isostructural conducting oxide SRO, might ex- 
plain the enhanced interfacial dielectric properties of the 
former electrode material, and might also help shed some 
light on the qualitatively different behavior of PTO/Pt 
and BTO/Pt. 

To investigate this hypothesis, we start by looking at 
the relaxed geometry of the interfacial Pt2 and AO lay- 
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FIG. 3: Schematic representation of the sahent features of the 
Pt2-AO interface. Only interfacial atoms are shown; empty, 
red, and blue circles denote Pt, O, and A atoms respectively, 
(a) Unsupported film and metal surfaces. A small surface 
rumpfing is present only in the AO surface layer while the 
metal is fiat, (b) Centrosymmetric interface geometry, with 
significant rumplings 5, especially in the Pt2 layer, (c) Pro- 
posed mechanism for interfacial ferroelectricity, driven by the 
unstable Pt-0 bond. Right: example of a two-body binding 
energy curve having negative curvature at d(0). 



ers when the film is in the centrosymmetric configura- 
tion, schematically shown in Fig. [31 The significant rum- 
pling 5 which is present in both Pt2 and AO (values are 
reported in Tab. |TT]) brings Pt(l) and the interface O 
atoms in much closer contact than the neighboring Pt(2) 
and A cation, indicating that only the former two atoms 
form a true chemical bond. Moreover, the strong buck- 
ling in the Pt2 layer is indicative of a frustrated bonding 
environment, where the repulsive A-Pt interaction is in 
competition with the attractive O-Pt bond. As a result, 
the O-Pt distance d in the relaxed interface structures is 
in both cases significantly larger than the typical Pt-0 
equilibrium distance of about 2.0 A, which has been de- 
termined experimentally and theoretically in a number of 
bulk oxide phases [2^. Comparing the interface systems, 
d is even larger in BTO/Pt than in PTO/Pt, indicating 
that the Pt-0 bond is weaker in the former system. 

To quantify this difference in strength, we examine 
the local force constants characterizing the Pt-0 inter- 
facial bond. At the PTO/Pt interface we have Kp^_Q = 
— 0.039 a. u., which is comparable to the typical inter- 
atomic force constants in oxide crystals and is consistent 
with a stable chemical bond. Strikingly, at the BTO/Pt 
interface we have Kp^._Q = 0.001 a.u., a very small and 
positive value, which indicates that the Pt-O bond here 
is slightly unstable (i.e. there is no restoring force asso- 
ciated to it). This indicates that the structural frustra- 
tion present at the interface is strong enough in BTO /Pt 
to pull the Pt-0 bond into an unstable regime [see Fig. 
3 (right)]. This qualitative difference between PTO/Pt 
and BTO /Pt is entirely responsible for their inequivalcnt 
polar response, as we shall see in the following. 

In order to link the stability of the Pt(l)-0 bond to 
the dielectric response of the capacitor, we perform a 
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FIG. 4: Ionic displacements induced by a small field of D = 
10"^a.u. (0.455 ^C/cm^). (a) PTO/Pt; (b) PTO/Pt modi- 
fied by artificial addition of force constant k = — 0.05a.u.; (c) 
BTO/Pt. The color/symbol scheme is the same as in Fig. [2] 
The horizontal axis is the x coordinate, normal to the inter- 
face plane, and refers to the relaxed atomic positions in the 
centrosymmetric geometry. 



computational experiment in which we artificially mod- 
ify the force constant matrices K-^ of PTO/Pt and 
BTO /Pt by adding a negative harmonic term of the form 
Eh = (fc/2)[a;o-a;pt-d(0)]2 between Pt(l) and O, which 
weakens the corresponding bond without changing the 
centrosymmetric equilibrium geometry. Our goal is to 
demonstrate that, by choosing an appropriate value of fc, 
one can destabilize the Pt(l)-0 bond at the PTO/Pt 
interface and therefore reproduce the behavior (dead- 
layer thickness and atomic displacement pattern) of the 
BTO/Pt system. To that end, we choose k = — 0.05a.u. 
(a value slightly larger than K^^ q) and we recalculate 
the effective dead-layer thickness for PTO/Pt by insert- 
ing this modified K-^ matrix into Eq. ([7]). Remarkably, 
due to the additional harmonic term Eh, Ndcad — 2.71 
jumps to a very small value iVdoad = 0.08, practically 
identical to that of the unmodified BTO/Pt case. 

To further prove that the "modified" PTO/Pt behaves 
in all important respects like BTO/Pt, we plot in Fig. [4] 
the induced ionic displacements in three different cases: 
PTO/Pt (fc = 0), PTO/Pt (fc = -0.05), and BTO/Pt 
{k = 0). Note that the ionic responses to a small D 
field are very different in the first and last cases. In 
this context, it is remarkable that the middle "modified 
PTO/Pt" case is strikingly similar to BTO/Pt, even in 
many fine details of the induced relaxations at the inter- 
face layers. Deeper into the oxide, the displacements re- 
flect the different ferroelectric mode patterns (e.g., there 
is a larger displacement for Pb than for Ba because of 
the stereochemical activity of the Pb lone pairs), but the 
behavior of the first oxide layers is nearly identical. This 
unambiguously demonstrates that the dissimilar behav- 
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FIG. 5: Computed local properties at the BTO/Au inter- 
face, (a) Ionic displacements induced by a small field of D = 
0.001 a.u.; (b) Inverse permittivity profile (thick solid black 
curve) compared with corresponding BTO/PT results (thin 
solid red curve) as in Fig. [T] The thin dashed lines in (b) 
indicate the values of the bulk inverse permittivity of BTO 
at D = for two different values of the in-plane strain, 0% 
(black) and —2.2% (black), corresponding to the cubic lattice 
constants of BTO and STO, respectively. 



ior in the BTO/Pt and PTO/Pt cases is almost entirely 
due to the change in stability of the Pt-0 bond. 

Now, to have a mechanism for ferroelectricity, a struc- 
tural instability needs to be coupled to the polarization. 
Comparing the original and the "modified" PTO/Pt 
ionic displacement patterns in Fig.[3{a) and (b), it is ap- 
parent that the stretching of the Pt-0 bond is associated 
with i) a shift of the whole ferroelectric slab in the oppo- 
site direction to D, and ii) with a combined rumpling of 
the AO and Pt2 layers. Since the atoms in both Pt2 and 
AO carry significant values of (see Fig. [2]), both effects 
have a strong impact on the polarization. We therefore 
identify this mechanism as a novel form of "interfacial 
ferroelectricity." Here we have a fiagrant breakdown of 
phenomenological models based on bulk properties of the 
parent materials. 

These arguments lead to a quite general prediction, 
namely that AO-terminated perovskite ferroelectrics can 
show a strong interfacial enhancement of the ferroelectric 
properties when weakly bonded to a simple metal. To 
test this prediction, we performed a supplementary cal- 
culation for a BTO-based capacitor by replacing the Pt 
electrodes with a more inert material, Au. In this case 
the equilibrium lattice constant of bulk BTO was used 
for the epitaxial strain constraint, to better match the 
experimental setup of Ref . [1^ , but otherwise all compu- 
tational parameters were kept the same. As before, we 
first relaxed the centrosymmetric geometry and then cal- 
culated the ionic and electronic response to a small field. 
The calculated induced ionic displacements for BTO/Au, 



6 



presented in Fig.[5lja), show the same quahtative features 
i) and ii) discussed above, but their magnitude is about 
four times larger than in BTO/Pt, reflecting the weaker 
Au-oxide interaction. Consistent with the enhanced dis- 
placements, we obtain a negative dead-layer thickness of 
-^dead — —5.5 layers, which manifests itself as an unusu- 
ally large dip in the local permittivity profile, shown in 
Fig. ISlh). (Note the rather dissimilar bulk "ferroelec- 
tric strength" of BTO at the different in-plane strain 
states.) In other words, the actual capacitor with 8.5 lay- 
ers of BTO and real Au electrodes has the same degree of 
ferroelectric instability as ~19.5 layers of BTO between 
"ideal" electrodes! It is quite remarkable that such a 
huge effect can be caused by only two atomic monolayers 
at the interface (Au2 and BaO). 

In conclusion, we have used a rigorous finite-field ap- 
proach to study the influence of the electrode-film inter- 
face on the ferroelectric instability of nanoscale BaTiOs 
and PbTiOa capacitors. At the AO-terminated interfaces 
with simple metals, we demonstrate a strong correlation 
of the overall ferroelectric response of the capacitor to the 
strength of the interfacial metal-oxygen bond. In cases 
where this bond is especially weak, such as the BaO-Au 
interface, we find an unusually large enhancement of the 
ferroelectric instability. This result is in striking con- 
trast with the conclusions of the Thomas-Fermi screen- 
ing model, demonstrating that a microscopic analysis is 
generally necessary to describe interfacial effects in fer- 
roelectric capacitors. This mechanism also suggests a 
possible route to the fabrication of nanoscale devices free 
from interfacial size effects. 



tries (a mirror symmetry plane is first imposed to relax 
the geometry in the paraelectric state, and later relaxed 
when calculating the linear response to a displacement 
field). 

The insulating character of the system along the stack- 
ing direction (the insulator is thick enough so that direct 
tunneling is suppressed) allows for a rigorous definition 
of the macroscopic polarization and its coupling to an 
external bias voltage 16, 19]. Because at the quantum- 
mechanical level all wavefunctions are mutually coupled 
and there is no clear way to separate free conduction 
charge from the bound polarization charge we adopt the 
convention of treating all charges as bound charges, so 
that £) is a constant throughout the heterostructure f27| 
just as for a purely insulating system 28|. What would 
customarily be called the macroscopic free charge cr on 
the top of the electrode is here reinterpreted as a bound 
charge associated with a polarization P = cr in the metal- 
lic region. (Indeed D = Aira, since the electric field must 
vanish in the metallic region.) 

All quantities requiring a derivative with respect to D 
were performed by finite differences between a calculation 
with mirror symmetry imposed {D — 0) and a fixed-D 
calculation with D = 0.001 a.u. The longitudinal force 
constant matrices and the dynamical charges were com- 
puted by finite differences by taking displacements of 0.4 
mA along the tetragonal axis. 
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Methods 

Our calculations are performed within the local- 
density approximation of density-functional theory and 
the projector-augmented- wave method |l1S], with a 
planewave basis cutoff energy of 40 Ry. We constrain 
the in-plane lattice constant to the calculated theoreti- 
cal equilibrium value for cubic SrTiOa (a = 7.276 a.u.) 
unless otherwise specified. We set the thickness of the 
insulating film to N unit cells (the actual values are re- 
ported in Table HI . The thickness of the metal electrode 
slab was 7.5 unit cells of SRO or 11 atomic layers of Pt, 
which were sufficient to converge the interface properties 
of interest. We use a 6 x 6 Monkhorst-Pack sampling 
of the surface Brillouin zone, which reduces to 6 spe- 
cial points within the tetragonal symmetry. A Gaussian 
smearing of 0.15 eV was used to accelerate convergence 
of the Brillouin-zone integrations. We carefully checked 
that the band alignment at the interface and the choice 
of the smearing scheme is such that any spurious pop- 
ulation of the conduction band of the insulating film is 
avoided. All the structural degrees of freedom were fully 
relaxed subject to the translational and point symme- 
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